Band structure calculations for Ba 6 Ge 2 5 and Ba 4 Na 2 Ge 2 5 clathrates 
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Electronic band structures for BaeGe25 and Ba4Na2Ge25 clathrates are calculated using linear 
muffin-tin orbital method within the local density approximation. It is found that barium states 
strongly contribute to the density of states at the Fermi level and thus can influence the transport 
properties of the compounds. A sharp peak of the density of states is found just at the Fermi level. 
It is also shown that the shifting of barium atoms toward experimentally deduced split positions in 
Ba6Ge25 produces a splitting of this peak which may be interpreted as a band Jahn- Teller effect. 
If the locking of the barium atoms at the observed structural phase transition is assumed, this 
reduction of the density of states at the Fermi level can account for the experimentally observed 
decrease of the magnetic susceptibility and electrical resistivity at the phase transition, and the 
values of density of states are in agreement with low temperature specific heat measurements and 
variation of superconducting transition temperature with pressure. 

PACS numbers: 71.20.-b,71.20.Tx 



I. INTRODUCTION 

BagGe25 clathrateBBi attracted considerable scientific 
interest and initiated an intense experimental activity be- 
cause of its rich physical properties. This class of com- 
pounds is investigated primarily because of their ther- 
moelectric properties!] The concept of PGEC, phonon- 
glass electron-crystal, was introduced. Massive alkali and 
alkali-earth metal atoms are expected to rattle in over- 
sized cages usually formed by Si, Ge or Sn atoms. Thus, 
they scatter the lattice phonons and reduce the thermal 
conductivity. On the other hand, electrical conductivity 
is believed to be more constrained to the cage and not 
much affected by rattling. 

BagGe25 undergoes a two-step first order structural 
phase transition at temperatures 215 and 180 K.ErEI The 
structural phase transition is accompanied by dramatic 
change of transport properties. Above the transition it 
shows a metal-like conductivity with fairly large value 
for this class of compounds. At the phase transition 
this behaviour is profoundly changed; resistivity abruptly 
increases, then slowly rises with lowering the tempera- 
ture and saturates at low temperatures. The Hall coef- 
ficient smoothly decreases with decreasing the tempera- 
ture and is not remarkably affected by the phase transi- 
tion. Magnetic susceptibility is negative and temperature 
independent above and abruptly decreas^.at the phase 
transition similarly to the thermopowerfl'Q It has been 
observed fiat hydrostatic pressure suppresses thenphase 
transition.0 The homologue clathrate Ba4Na2Ge25 ,B with 
smaller lattice parameter, does not undergo a phase tran- 
sition and shows typical metal-like behaviour.El At low 
temperatures bpth clathrates experience superconduct- 
ing transitions. 

In this article we consider the electronic band struc- 
ture of the compounds and the influence of the structural 
phase transition on the electronic band structure as well 
as the possible consequences on the transport properties. 
In Section II some structural details are considered. De- 



tails of the calculations are described in Section III. The 
results and the discussion are presented in Section IV. 
Finally, the summary and the conclusions are given in 
Section V. 



II. STRUCTURE 

The structure of BagGe25 belongs to the space group 
P4i32 (No. 213). There are 124 atoms in the unit cell. 
All the structural details can be found in Ref. |T|,[s|,|9||lC|. 
Thorough description of the structure can also be found 
in Ref. [ll],[L2] on structurally similar Sn clathrates. Ger- 
manium atoms form the framework which defines three 
different types of cages filled with barium atoms. The 
main building blocks are distorted germanium pentagon- 
dodecahedra with Bal atoms inside. They form a chiral 
network with large cavities, where Ba2 and Ba3 atoms 
are placed. In the Ba4Na2Ge25, two sodium atoms are 
randomly distributed in three Ba2 sites per formula unit. 
In total, 17 of 25 germanium atoms per formula unit are 
four- fold coordinated by other germanium atoms. The 
remaining 8, usually denoted as Ge3 and Ge5, are three- 
fold coordinated. The average distance between germa- 
nium atoms (~2.54 A) is somewhat larger than for Ge 
in diamond type structure (~2.45 A) and the angles are 
distorted. Typical nearest distance between barium and 
germanium atoms is around 3.5 A. The nearest distance 
between barium atoms is around 4.70 A for Ba2-Ba3 in 
undistorted structure (see also Ref. |l0| ). 

Due to highly aspherical electron density deduced from 
X-ray diffraction data, split positions for Ba2 and Ba3 
atoms are assumed. The most pronounced effect at the 
structural phase transition is the increase of the atomic 
displacement-parameter for Ba2, and somewhat less for 
Ba3 atom.Bii3 It is believed that below the phase tran- 
sition Ba2 and Ba3 atoms are actually locked in the 
split sites. However, since no superstructure is observed, 
they should be randomly distributed among split sites. 
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The germanium cage also undergoes modifications at 
the phase transition and there is not only a difference 
between the room temperature and the low tempera- 
ture structure, but also the variation pwith temperature 
through the whole temperature range. t3 



III. COMPUTATIONAL DETAILS 

Band structure calculations were performed within Lo- 
cal Density Approximation (LDA) using Linear Muffin 
Tin Orbital (LMTO) method in Atomic Sphere Approxi- 
mation (ASA) with combined correction terms taken into 
account oll3 Scalar relativistic version was used with von 
Barth-Hedin exchange-correlation potential.ta 

The calculations for the room- and the low- 
temperature structural data for Ba 6 Ge25 were performed. 
The atomic parameters for different models used in cal- 
culations are given in Table |. The barium atoms were 
placed in their average positions within higher symmetry 
model, i.e., space group P4i32. In Ba4Na2Ge25, the ran- 
dom distribution of the Na and Ba atoms at Ba2 site was 
neglected and the calculations were performed for lower 
symmetry model, namely, group P4i2i2 (No. 92), with 
sodium atoms at one of the 86 symmetry sites. 

The calculations with Ba2 and Na atoms shifted from 
their average positions in the original space groups 
(No. 213 and No. 92, respectively) were performed using 
lower symmetry models in space groups P2 1 3 (No. 198) 
and P4i (No. 76), respectively. In this way the assumed 
randomness of the locking was not considered. 

Most of the results presented in this work were ob- 
tained from self-consistent calculations on 12 x 12 x 12 
k-space mesh. For the detailed band structure around 
Fermi level the calculations were performed for the 16 x 
16 x 16 mesh. The values of the DOS at the Fermi level 
for these cases showed corresponding variations of about 
1 eV^ 1 per unit cell. The total energy was negligibly 
affected by this change of the mesh size. 

The atomic sphere radii were determined so that their 
overlap was minimised. This resulted in the radii of 
about 1.5 A (2.9 a.u.) for germanium and 2.5 A (4.8 a.u.) 
for barium atoms. It is worth to note that, as the radii of 
germanium spheres are rather small, germanium p-states 
extend far beyond the sphere boundaries. The large part 
of these states, so called "tails," reside in the neighbour- 
ing barium spheres. This should be taken into account 
when considering barium and germanium charges which 
are defined as the integrals of the electron density over 
the corresponding sphere. They are clearly identified in 
the empty spheres of hypothetical n6Ge25, where the bar- 
ium atoms are replaced by the empty spheres of the same 
radii. In this case a high electron density near the bound- 
ary can be observed which diminishes toward the center 
of the empty sphere. Because of the large barium spheres, 
also the states of /-symmetry were used in the expansion 
of the wave function. In the germanium spheres the wave 
function was expanded up to the states of <i-symmetry. 
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FIG. 1: The electronic band structure and the density of 
states for Ba6Ge25 in the low temperature structure at 10 K 
with Ba atoms in the average positions within space group 
P4i32. 



IV. RESULTS AND DISCUSSION 

Here we present and discuss the results for the low 
temperature structure of Ba 6 Ge25 at 10 K, first with Ba 
atoms in the average sites. We focused the analysis on the 
density of states and the bonding properties between the 
barium atoms and the germanium cage. No distinctive 
qualitative difference in the DOS between the room and 
the low temperature structures within higher symmetry 
model was found. Further, we investigated the effects of 
shifting the barium atoms toward split positions on the 
band structure. This was done within the low tempera- 
ture structure, because there the split positions are the 
most distant and the locking of the barium atoms is as- 
sumed. Finally, the results for the Ba4Na2Ge25, which 
does not undergo phase transition, are discussed. In par- 
ticular, the influence of the sodium on the band structure 
is examined in comparison with barium. 



A. Electronic band structure and density of states 

The electronic band structure with its density of states 
calculated for BaeGe25 in the low temperature undis- 
torted structure are shown in Fig. [fl The overall pic- 
ture is similar for the room temperature structure. Due 
to the large number of atoms in the unit cell there is 
a large number of bands and a complex electronic band 
structure. 
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TABLE I: Wyckoff positions splitting from the group No. 213 in its subgroups. The average structure of BaeGe25 belongs 
to the space group P4i32 (No. 213). The values for the atomic coordinates were taken from Ref. [l]Jto|. Shifting the Ba2 
atom toward one of the split positions of Ba4 site demands lowering the symmetry to the space group P2i3 (No. 198). For 
Ba4Na2Ge25 (atomic coordinates were taken from Ref. ^) the symmetry had to be lowered to the group P4i2i2 (No. 92) so 
that Na atoms could be orderly placed on Ba2 sites. Shifting the Ba or Na from Ba2 site demands lowering the symmetry to 
the space group P4i (No. 76). 



Sites 


No. 213 


No. 198 


No. 92 


No. 76 


Gel 
Ge5 
Ge6 
Ba4 (Na) a 


24e 


x, y, z 


126 


x, y, z 


86 


x, y, z 


4a 


x, y, z 


4a 


-x+1/2, y+1/2, -z+1/4 


86 


z+7/8,x+3/4,y+3/8 


4a 


z+7/8,x+3/4,y+3/8 


126 


y+3/4, x+1/4, -z+1/4 


4a 


-z-3/8, x+1/4, -y-1/8 


86 


y+1/4, z+5/8, x+1/8 


4a 


y+1/4, z+5/8, x+1/8 


4a 


-y+1/4, z+1/8, -x+1/8 


Ba2 (Na) 
Ge4 


12d 


1/8, y, y+1/4 


126 


1/8, y, y+1/4 


86 


7/8, y+1/2, y+1/8 


4a 


7/8, y+1/2, y+1/8 


4a 


-3/8, y, -y+1/8 


4a 


y+3/4, y+3/4, 


4a 


y+3/4, y+3/4, 


Bal 

Ge2, Ge3 


86 


X, X, X 


4a 


X, X, X 


86 


x+3/4, x+1/2, x+7/8 


4a 


x+3/4, x+1/2, x+7/8 


4a 


-x+3/4, -x+3/4, -x+3/4 


4a 


-x-1/4, x, -x-5/8 


Ba3 


4a 


3/8, 3/8, 3/8 


4a 


3/8, 3/8, 3/8 


4a 


1/8, 7/8, 1/4 


4a 


1/8, 7/8, 1/4 



"Split positions of the Ba2 and Na. 



In Fig. |2|the /-projections of DOS in germanium atomic 
spheres are shown: Ge in diamond structure, hypothet- 
ical Ge-graphite, □eGe25, and Ba 6 Ge25 compounds are 
compared. In the DOS, s-, sp-, and p-bonding regions 
can be recognised in filled and empty clathrate similar 
to those in Ge-diamond or graphite. This similarity is 
to be expected since locally 17 of 25 germanium atoms 
per formula unit are four-fold coordinated (4b) and form 
covalent sp 3 -like bonds, like in the diamond structure. 
The other 8 germanium atoms are three-fold bonded (3b) 
more like in graphite structure. The total width of the 
bonding region is about 0.9 eV narrower in the empty and 
filled clathrates. The fundamental gap between bonding 
and anti-bonding region is larger in n6Ge25 than in the 
diamond structure by the same amount of 0.9 eV and 
there appears a gap between sp- and p-like bands in the 
valence region. These features are common for the other 
Ge and Si clathrates (e.g., Ref. The appearance of 
the gap in the valence region by its contraction was asso- 
ciated with the 5pr.kig patterns (pentagons) of germanium 
or silicon-atoms,t3't£l but there was also some criticism of 
the idea.El According to the latter, the small angular dis- 
tortion of the tetrahedrally bonded framework may also 
play an important role. While for Ge clathrates of type 
I and II an overlap of the s- and sp-region is observed,E3 
here one can see a clear separation. One should keep 
in mind that in BagGe25 structure, in contrast to the 
clathrates I and II, there are three bonded germanium 
atoms. However, the band structure calculation on hy- 
pothetical Ge-graphite, where all germanium atoms are 
three-fold coordinated, did not reveal pronounced ten- 
dency toward such a separating of the s- and sp-regions, 
as seen in Fig. [| 

The main difference in the site projected DOS for ger- 
manium atomic spheres between the empty and the filled 
clathrate is the reduction of the gap between bonding and 
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FIG. 2: Projected density of states from Ge atomic spheres in 
Ge-diamond, hypothetical graphite structures, BaeGe25 and 
DsGe25. For the clathrate compounds, the average contribu- 
tion from all germanium is shown. 



anti-bonding region in the filled clathrate. This immedi- 
ately shows the importance of barium atoms not only as 
the pure electron donors to the conduction band formed 
of germanium anti-bonding states. 

The three-bonded germanium atoms, six Ge5 and two 
Ge3, form a distorted cube around Ba3. Ba2 atoms lie on 
the lines passing from the central Ba3 through the cube 
faces formed by three Ge5 and one Ge3 atom. In this way 
they form a large octahedron around Ba3. The closest 
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TABLE II: Total charges in the atomic spheres and their radii 
for BagGe25 an d O e Ge25- 



At. Sphere 


Ba B Ge 25 


□ 6 Ge 25 


difference 


AS radius (a.u.) 




Bal, Dl 


4.613 


2.853 


1.760 


4.7175 


0.4 


Ba2, D2 


3.474 


1.864 


1.610 


4.6919 




Ba3, D3 


3.711 


1.776 


1.935 


4.8720 


0.2 


Gel 


3.520 


3.462 


0.058 


2.8663 


Ge2 


3.821 


3.805 


0.016 


2.9607 


0.0 


Ge3(3b) 


3.536 


3.403 


0.133 


2.9829 


Ge4 


3.613 


3.619 


-0.006 


2.8728 




Ge5(3b) 


3.436 


3.315 


0.121 


2.8675 




Ge6 


3.559 


3.499 


0.060 


2.8648 


0.1 



distances are between Ba2 and Ge3, 3.35 A, and Ba3 and 
Ge5 3.38 A. Ge3 atoms are somewhat more distant from 
Ba3, at 3.89 A. This environment and the nature of the 
bonds between Ba3, Ge3, Ge5, and Ba2 should be the key 
to understand the physical properties of the compound 
and its phase transition.!!^ 

Let us consider the charges in the atomic spheres. As 
mentioned in Section III, there is a large contribution of 
germanium states in the oversized barium spheres. This 
can be seen in Table || in the column for the charges 
of empty clathrate, n6Ge25. The large charges in the 
empty spheres are coming from the tails of germanium 
states. The highest value is in Dl. This is expected be- 
cause Dl has the closest Ge neighbours forming distorted 
pentagon-dodecahedron. Despite of the fact that D3 has 
the largest sphere, relatively smaller charge is found in- 
side due to the fewer Ge neighbours. For three-bonded 
germanium atoms it can also be noted that they have 
somewhat smaller charges inside their spheres in compar- 
ison to the four-bonded germanium, although e.g., Ge3 
has the largest atomic sphere from all germanium atoms. 
Among germanium atoms, three-bonded Ge3 and Ge5 
have the largest increase of the charge upon inclusion of 
barium atoms. This is in accord with the expectation 
of Zintl concept, namely, that three-bonded germanium 
atoms accept an electron from barium atoms to form the 
lone electron pairs. However, because of the small germa- 
nium spheres the absolute value of the difference is not 
so pronounced. Among barium spheres, the largest in- 
crease of the charge, in comparison to the corresponding 
empty spheres, undergoes the Ba3 sphere. We interpret 
this as the influence of the lone pairs, reflected in the 
Ba3 sphere as the tails from neighbouring three-bonded 
germanium atoms. This is illustrated in Fig. ^| where 
the differences of the ^-projected radial charge densities 
in the several atomic spheres are shown. The large elec- 
tron density near the boundary of Ba3 atomic sphere may 
be associated with the lone electron pairs of the neigh- 
bouring three-bonded germanium atoms. It was also 
shown by electron localisation function (ELF) aaalysis 
that these lone pairs are situated around Ba3 siteJij The 
oscillating part of the density, i.e., its difference shown 
in Fig. |^, may be interpreted as more atomic-like barium 
states. The biggest contribution comes from the states 




r(a.u.) 



r(a.u.) 



FIG. 3: Electron density difference between filled and empty 
clathrate in particular atomic spheres. 



of d-symmetry. 

We conclude that there is only a small charge transfer 
from barium to the cage, just in agreement with conclu- 
sion made for BasGai6Ge3o in Ref. ^l], with difference, 
that here all the arguments point out that a consider- 
able amount of the electron charge remains in the states 
with significant barium character. This is in contrast to 
the situation reported for K 8 Ge4 6 in Ref. ^2] where the 
authors conclude complete charge transfer from potas- 
sium atom to the frame. Both of these cited results were 
critically re-considered and charge transfer redefined in 
the recent paper .E3 However, as may be seen in Fig. |^, 
the large amount of the charge in the barium spheres re- 
sides near the boundaries, so that it may be said that the 
conduction electrons are more constrained to the cage. 

In the band structure calculations for the similar 
KeSn25 and KgSn23Bi2 compounds, a lifting of the three 
bands, originated only from three-bonded tin atoms, 
from the top of the valence region was found and ascribed 
to the interaction,-of the lone pairs formed by three- 
bonded tin atoms.Ej. But, as was pointed out there, in 
BaeGe2iIn4 this characteristic is absent due to a smaller 
spatial extension of germanium orbitals. We did not no- 
tice significant difference of DOS between three-bonded 
(Ge3, Ge5) and four-bonded germanium atomic spheres. 
However, as discussed above, we expect that the lone 
electron pair states mostly reside as tails in D3 atomic 
sphere. Therefore one might expect their contribution 
to the DOS in D3 sphere in n6Ge25. Indeed, we find a 
peak in DOS on the top of the unoccupied valence re- 
gion (Fig. ^) appearing only in D3 sphere. This peak 
remains in the filled clathrate, i.e., in Ba3 sphere, what 
can also be seen in Fig. |I| For Ba2 and Ba3 spheres a 



5 



Ba2 in Ba 6 Ge 25 
Ba2 s 

Ba2 p 1 

Ba2 d " 




Ba3 in Ba 6 Ge 25 

Ba3 s 

Ba3 p 

Ba3d 


J 

<•' ■■-..-*--,'. 


□ 2inD 6 Ge 25 

- 2s 

U2 p 

02d 




□ 3 in CeGegg 

- 3s 

3 p 

D3d 



-6 -4 -2 2 4 6 8 10 1B -4 -2 2 4 6 8 10 



Energy (eV) Energy (eV) 

FIG. 4: The projected density of states from atomic spheres 
of Ba2 and Ba3 in BaeGe25 and the equivalent empty spheres 
of D2 and 03 in □ 6 Ge 25 . 



large increase of DOS, compared to the empty D2 and D3 
spheres, in the region around the Fermi energy is found. 
The number of states is increased and this suggests that 
the barium states participate in the conduction band for- 
mation. 

This iain contrast to the results for the BagGaieGeso 
clathrateO where the inclusion of barium atoms inside 
the Ga-Ge cage left the band structure almost intact- 
This is also different to the cases of some Si clathratcscfl 
where the main effect of the inclusion of alkali metal was 
shifting of the particular state at the bottom of the con- 
duction band, but still predominantly of the cage-like 
character so that it was possible to neglect the overlap of 
the orbitals of guest atoms with host-cage. Here, how- 
ever, the distances between the barium and germanium 
atoms are shorter and stronger hybridisation may be ex- 
pected. 

Another important feature of DOS is that it shows a 
peak just at the Fermi level (Fig. Its value is around 
50 states/eV per unit cell in the low temperature struc- 
ture. The largest contributions come from barium atomic 
spheres mainly with e?-symmetry and from p-states in 
germanium atomic spheres. The shape of the DOS com- 
ing from barium and germanium spheres follows closely 
each other. This indicates a strong hybridisation of bar- 
ium and germanium states. The contribution from Bal 
spheres is significantly lower than that from Ba2 and 
Ba3 spheres, although Bal has the nearest Ge neigh- 
bours. This complies with the assumption that the Ba2 
and Ba3 atoms with the surrounding three-bonded ger- 
manium atoms determine the most unusual properties of 
the system. It is quite interesting that the similar peak 



of the DQS-at the Fermi energy was found for Ba§Si46 
clathrateE§E3 as well. In that case also the relatively 
shorter distances between the barium and silicon atoms 
were found (~3.3-3.5-A). We have also performed cal- 
culation for BagSi25.c3 The lattice parameter is smaller 
compared to the Ge clathrate and the fundamental gap 
disappears completely but the peak in the DOS at the 
Fermi level with strong contribution of the states with 
d-symmetry coming from the barium atomic spheres re- 
mains (Fig. ||). It is also interesting that the subband- 
gap between s- and sp- regions (not shown in the Fig. |5|) 
appears in the BagSi25 as well. 

For a better understanding of the above results, i.e., 
the role that is played by barium states, we have per- 
formed the calculation with barium atoms only, BagC^s, 
without the Ge framework. The distances between Ba2 
and Ba3 atoms of about 4.70 A are comparable to the dis- 
tances between atoms in Ba bcc structure with 4.34 A. 
The calculation showed the dominant contribution of d- 
symmetry states at the Fermi level. The s-states are split 
into bonding and the anti-bonding regions. The bonding 
states lie below Fermi level and anti-bonding above, so 
that they do not give significant contribution to the DOS 
at the Fermi level. The tails of the barium d-states in the 
empty germanium spheres are of s- and p-symmetry, so 
that the appreciable overlap of these states with germa- 
nium s- and p-states in the filled clathrate might be ex- 
pected. Concerning the relative energy parameters of the 
barium and germanium muffin-tin orbitals, they showed 
qualitative agreement with the values used for the atomic 
orbitals in LCAO calculations .E3 All these facts suggest 
that the appearance of the barium <i-states and their hy- 
bridisation with germanium p-states at the Fermi level is 
quite plausible. 

The possible importance of the hybridisation of the 
Ba 5d states with Si cage was also pointed out for 
BasSi46.li3 The results of the experiment with ultrahigh- 
resolution photoemission spectroscopy for BasSi46 sug- 
gest even stronger hybridisatkp of Ba and Si states than 
obtained in LDA calculation.E3 NMR experiments on this 
compound could be interpreted if strong contribution p£ 
Ba (i-states to the DOS at the Fermi energy is assumed.ESl 



B. Band structure near the Fermi level versus Ba2 

shift 

The presence of the peak in the DOS at the Fermi 
level suggests a possible instability of the system. Since 
there are strong indications of locking of Ba2 atoms at 
the phase transition, as mentioned in Section II, we in- 
vestigated the modifications of the band structure due 
to the shifting of Ba2 atoms within the low temperature 
structure. 

The most interesting effect of this shifting of Ba2 
atoms, performed as explained in Section III, is the split- 
ting of the peak at the Fermi level and, consequently, the 
reduction of DOS. The effect was found to be the most 
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FIG. 5: Contribution to DOS from all Ba d-states, all Ge 
p-states and all three-bonded Ge p-states per formula unit. 
The total DOS is also shown. In the figure below the same is 
shown for the case of BagSi25, with Si atoms instead of Ge. 



pronounced with Ba2 at 80% distance from the higher 
symmetry position toward the experimental split posi- 
tion. In Fig. ||, the band structure together with the DOS 
in the close vicinity of the Fermi level for this structure 
is compared with the higher-symmetry (average) struc- 
ture. The effect may be considered as a band Jahn- Teller 
effect, though in the band structure along the symmetry 
directions no obvious removing of the degeneracy is ob- 
served. The equivalent scenario was proposed for some 
other compounds with somewhat similar metal-insulator 
transitionO The conventional Jahn- T oller eff ect was con- 
sidered also for some other clathrates.BilESEil 

The total energy was also compared for shifted posi- 
tions. The minimum was found to be at the off aver- 
age position. The potential energy surface for Ba atoms 
showed binding stronger than expected, but the variation 
of the total energy with respect to the shifting of atoms 
is known not to be highly accurate within the atomic 
sphere approximation. That is also the reason why the 
structure optimization was not performed. 

One has to be careful with the interpretations, espe- 
cially because the whole framework is slightly changing 
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FIG. 6: Electronic band structure and the density of states 
around the Fermi energy for the Ba2 at the average position 
and at 80% toward the split site in BasGe25. Dashed lines 
designate the bands that are crossing the Fermi level. 



with the temperatuKe-, and is strongly disturbed during 
the phase transition.liJ As mentioned before, we did not 
notice a qualitative difference in the DOS between the 
room temperature and the low temperature higher sym- 
metry structures. For both of them the peak at the Fermi 
level is found. The splitting of the peak when Ba2 atom 
is shifted toward one of the split positions indicates that 
the locking of barium atoms might be the most important 
process at the phase transition. If we assume that at the 
phase transition this locking of barium atoms occurs then 
the reduction of the DOS at the Fermi level is in a quali- 
tative agreement with the observed decrease of magnetic 
susceptibility. The decrease of the value of DOS is found 
to be around 25 eV~ x per unit cell which corresponds to 
a drop of the Pauli susceptibility of 2.08 x 10 -4 cm^/mol 
per formula unit J&fA, agrees in order of magnitude with 
the observed one£li3 Moreover, the area of the Fermi sur- 
face is reduced because there are less bands crossing the 
Fermi level, what can also be seen in Fig. ^. As a con- 
sequence, electrical conductivity is expected to be lower 
and this is also in accord with the experimental obser- 
vation. It should be emphasised that we did not take 
into account the effect of the disorder of locked barium 
atoms. However, due to a very short electron mean free 
path, which is much smaller than the unit cell parameter, 
as estimated from experiments and a simple free electron 
modeljQ the barium disorder effects should not reflect 
themselves in the transport properties. The same ar- 
gument may be applied to understand experimental Hall 
coefficient results, i.e., the electrons might not be able 
to close the trajectory on the Fermi surface and there- 
fore the Hall coefficient might not be strongly influenced 
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by the phase transition. The absolute value of DOS at 
the Fermi level is in agreement with the low temperature 
specific heat measurements B If we assume the locking of 
barium atoms at low temperatures and take the corre- 
sponding value of DOS of 25 eV -1 at the Fermi level, 
we obtain 7 = 15 mJ/molK 2 per formula unit. The ex- 
perimental value is 21 mJ/molK 2 per formula unit. For 
the Ba4Na2Ge25, discussed in the next section, we obtain 
29 mJ/molK 2 per formula unit in comparison with the 
experimental 33 mJ/molK 2 per formula unit. The mea- 
surements of the variation of superconducting transition 
temperature under pressure show that T c is increased 16 
times, from the value of 0.24 K at zero pressure, as the 
structural phase transition is suppressed. Introducing the 
values of DOS for higher and lowes symmetry structure 
in the McMillan expression for Tji3 provides the ratio of 
21, in a relatively good agreement with the observed one. 



C. The results for Ba4Na2Ge25 

The total band structure for Ba4Na2Ge25 is qualita- 
tively similar to that found for BagGe25. The Fermi level 
is shifted downwards because there are less electrons and 
the gap between bonding and anti-bonding region is nar- 
rower since the lattice parameter is smaller compared to 
BaeGe25. The peak of DOS at the Fermi level is also 
found. 

In order to examine the influence of the sodium on the 
band structure, we performed a calculation on a fictitious 
Ba4.n2Ge25 compound and compare it with Ba4.Na2Gc25. 
The inclusion of sodium atoms introduces a lowering of 
the bands at the beginning of the conduction band to- 
gether with the most upper parts of the valence bands, 
so that the fundamental gap has nearly the same value. 
The conduction bands are modified. Although some sim- 
ilarities may be observed, rigid-band model could hardly 
be applied. In the valence region, however, apart from 
the bands just below the fundamental gap, near F point, 
the band structures are almost identical. Introduction of 
sodium does not influence the bonding states of germa- 
nium. The /-projected DOS in sodium atomic spheres 
are shown in Fig. |^, and compared with the pertinent 
empty sphere. In the sodium sphere a peak of mainly 
s-symmetry at the Fermi level is found, where the mod- 
ifications of the band structure, considered above, is ob- 
served. 

We have also performed a calculation on Ba4Ba2Ge25, 
where we have inserted barium instead of sodium atoms 
and kept the same symmetry and the structure. It is 
interesting that no shifting of the conduction band is ob- 
served. However, the band structure undergoes stronger 
changes than in the case of sodium. Even the valence 
region, i.e., bonding germanium states, is significantly 
affected. The /-projected DOS in the barium spheres, as 
in BagGe25, discovers the appearance of the states of d- 
symmetry (Fig. ||). The large contribution of (i-symmetry 
states are present around and above the Fermi level, in 
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FIG. 7: Comparison of the band structures around the 
Fermi level of the fictitious Ba4d2Ge25, Ba4Na2Ge25, 
and Ba4Ba2Ge25 with same structure and symmetry as 
Ba4Na2Ge25- 



contrast to the sodium and the empty spheres, where 
their contribution is small. 

The electron charges, projected on the states of par- 
ticular /-symmetry, in the sodium, barium, and empty 
sphere s fo r corresponding compounds are given in the 
Since in the empty spheres only the tails 



Table III 



from germanium states are present, the difference to the 
spheres of sodium and barium should give the contribu- 
tion of more atomic-like states. In the sodium sphere, 
the charge is mainly in the states of s-symmetry, but 
comparable contributions are found in the states of p- 
and somewhat less in cZ-symmetry states. Due to the 
change of the potential in the sodium sphere in com- 
parison with empty one, it is hard to disentangle the 
contribution coming from germanium tails. The values, 
however, indicate hybridisation of sodium s- and p-states 
with germanium states, although not so strong as in the 
case of barium in Ba4Ba2Ge25, where a distinctively large 
contribution of (i-symmetry states in the barium sphere 
is found, which again points to the strong hybridisation 
of barium cZ-states with germanium p-states. 

A shifting of sodium atoms to one of the split positions 
did not produce such a splitting of the band structure as 
in the case of BagGe25, although some reduction of DOS 
is observed. 
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TABLE III: The electron charge in the particular Ba, Na, and □ atomic sphere of Ba4Na2Ge25, Ba4D2Ge25, and Ba4Ba2Ge25, 
decomposed in the states of particular orbital angular momentum quantum number. The calculations for all compounds were 
performed within the same space group (No. 92) with the same structure parameters. Ba2, Na, and □ atomic spheres are all 
of the same radius of 4.59 a.u. and are equivalent to the Ba2 site in the Ba 6 Ge25 in the original space group No. 213. 



At. sphere 


Ba4Na2Ge25 
s p d f 


Ba4d2Ge25 
s p d f 


Difference 
s p d f 


Ba2 

Na,n 


0.53 0.76 1.52 0.43 
0.65 0.89 0.71 0.39 


0.49 0.76 1.51 0.43 
0.34 0.65 0.58 0.37 


0.04 0.00 0.01 0.00 
0.31 0.24 0.13 0.02 




Ba4Ba2Ge25 


Ba4d2Ge25 




Ba2,D a 


0.55 0.76 1.53 0.43 


0.34 0.65 0.58 0.37 


0.21 0.11 0.95 0.06 



"These correspond to the Na site in Ba4Na2Ge25 
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SUMMARY AND CONCLUSION 



Energy (eV) 

FIG. 8: The /-projected density of states from Ba2 and Na 
atomic spheres in Ba4Na2Ge25 and □ sphere, equivalent to 
the Na sphere, in Ba4Q2Ge25. They all correspond to the 
Ba2 site in Ba6Ge25 in the original space group No. 213. 



We have performed LDA band structure calculations 
for the two clathrates, BaeGe25 and Ba4Na2Ge25, using 
LMTO method within the atomic sphere approximation. 
The overall band structure reveals similarities with Ge in 
diamond type structure and the comparison with empty 
Ge cage shows very similar structure in the bonding re- 
gion. However, the conduction bands, in particular at 
the Fermi level, have a strong contribution of barium 
states. The analysis of the /-projected DOS shows that 
these states are of strong d-symmetry character. The cal- 
culations yield the peak of the DOS at the Fermi level 
that hints on the structural instabilities of the system. In 
BagGe25, shifting of Ba2 atoms, which have the largest 
and the most aspherical atomic displacement parameter, 
toward one of the experimentally deduced split positions 
splits the peak and reduces the DOS. The splitting of 
the peak at the Fermi level, due to the shifting of Ba2 
atoms and corresponding symmetry lowering, may there- 
fore be considered as the band Jahn- Teller effect. Barium 
in BaeGe25 clathrate can not be considered as the elec- 
tron donor only and this compound does not fit into the 
simple concept of PGEC. In the case of Ba4Na2Ge 2 5, the 
shifting of sodium atoms to one of the split positions does 
not produce such a pronounced tendency toward splitting 
the peak of DOS. If the locking of Ba2 atoms in BagGe25 
at the phase transition is assumed, the reduction of DOS 
at the Fermi level is in the agreement with experimentally 
observed abrupt decrease of the magnetic susceptibility 
and with the observed variation of superconducting tran- 
sition temperature with pressure. It may qualitatively 
account for the observed reduction of the electrical con- 
ductivity. The specific heat coefficient 7 obtained from 
the low temperature specific heat measurements is also 
in agreement with the calculated value of DOS at the 
Fermi for Ba@Ge25 with barium atoms locked, and for 
Ba 4 Na 2 Ge25 m the average structure. 
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